BRL  R 1749 


ALGEBRAIC  MANIPULATION 


M.  A.  Hirscnberg 


November  1974 


inoMwad  far  auhlle  mUiiii  dlitHhutlon  unllnrittd. 

• rr • ****  ■*’  r — ■ • - - 


USA  BALLISTIC  RESEARCH  LABORATORIES 

ABiRDiiN  MOVING  GROUND;  MARYLAND 


Destroy  this  report  when  it  is  no  longer  needed. 
Do  not  return  it  to  the  originator. 


Secondary  distribution  of  this  report  by  originating 
or  sponsoring  activity  is  prohibited. 

Additional  copies  of  this  report  nay  be  obtained 
from  the  National  Technical  Infoxmation  Service, 

U.S.  Department  of  Commerce,  Springfield,  Virginia 
22151. 


The  findings  in  this  reoort  are  not  to  be  construed 
an  official  Department  of  the  Army  position,  unless 
so  designated  by  other  authorized  documents. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


1.  REPORT  NUMBER 


BRL  REPORT  NO.  1749 


«.  TITLE  (end  Subtitle) 


READ  INSTRUCTIONS 

QOTAot?  PAlini  n'rrv r>  E>ntni 
u^i1  unit  cumruii  a in  lj  r vkwi 


3.  RECIPIENT’S  CATALOG  NUMBER 


S.  TYPE  OF  REPORT  ft  PERIOD  COVERED 


SYMBQLANG  - A SLIP  EXTENSION  FOR  ALGEBRAIC 
MANIPULATION 


6.  PERFORMING  ORG.  REPORT  NUMBER 


■ AUTHORfaJ 


8.  CONTRACT  OR  GRANT  NUMBERfaJ 


M.  A.  Hirschberg 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

USA  Ballistic  Research  Laboratories 
Aberdeen  Proving  Ground,  Maryland  21005 


10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  A WORK  UNIT  NUMBERS 

Program  Element  62618A 
Project  1T662618AH80 


11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  12.  REPORT  DATE 

US  Army  Materiel  Command  NOVEMBER  1974 

5001  Eisenhower  Avenue  is.  number  of  pages 

Alexandria,  VA  22304  29 


14.  MONITORING  AGENCY  NAME  a ADDRESSff/  different  inn  Controlling  Olllce)  IS.  SECURITY  CLASS,  (ol  thle  report) 

UNCLASSIFIED 

15a,  DECL.A.SSI  FI  cation/ DOWN  GRADING 
SCHEDULE 


17.  DISTRIBUTION  STATEMENT  (ol  the  abstract  entered  In  Block  20,  II  dlllerent  Irom  Report) 


19.  KEY  WORDS  (Continue  on  rover**  elde  //nec****ry  and  Identify  by  block  number) 


Algebraic  Manipulation 
Symbol  Manipulation 
List  Processing 
List  Manipulation 


Formula  Manipulation 
Data  Manipulation 
Formula  Translation 
Symbolic  Processing 
String  Processing 


20.  ABSTRACT  (Continue  on  rmveree  elite  If  ne  creme  ary  and  Identify  by  block  number) 


SYMBOLANG,  originally  developed  by  A.  Lapidus,  M.  Goldstein,  S.  Hoffberg,  and 
expanded  by  H.  Bernstein,  S.  Greenspan,  A.  Magnus,  and  others  is  a collection 
of  FORTRAN-cal lable  subroutines  which  can  perform  arithmetic  operations 


Ulhf  TSrf  1 Hill  at»  /I’Jiri  ci  nn  of  p 

juuli  aviiun  j iaai  tlpiiwuLiL/ii)  ur  v j.  oiuii  } w tv  • j • 


evaluations,  and  differentiations  on  expressions  represented  as  SLIP  lists. 
SYMBOLANG  consists  of  seventy-two  subprograms  and  is  operational  on  the 
BRLESC  II  computer. 


dd 


FORM 
J AM  73 


EDITION  OF  f NOV  6S  IS  OBSOLETE 


Unclassified 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Date  Entered) 


TABLE  OF  CONTENTS 


I  INTRODUCTION 


5 


A.  Background 5 

B.  Ordering  and  Simplification  7 

C.  System  Usage 7 


II  SYMBQIANG  SUBROUTINES 


A.  Basic  Set 9 

B.  Expanded  Capability  Set 21 

III  PROGRAM  EXAMPLES 22 


A.  Inputs  and  Outputs 22 

B.  Arithmetic  Calculation,  23 

C.  Differentiation 24 

D.  Truncation 25 

E.  Evaluation 26 

IV  CONCLUDING  REMARKS 28 

DISTRIBUTION  LIST 29 


3 


I.  INTRODUCTION 


A.  Background 


SYMBOLANG  was  written  by  Lapidus,*  Goldstein  and  Hoffberg  between 
1965  and  1967  and  extensively  expanded  by  Bernstein,2  Greenspan , Magnus 
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It  is  a collection  of  FORTRAN- call able  subroutines 
which  can  perform  arithmetic  operations  (addition,  multiplication, 
subtraction,  division,  etc.),  substitutions,  evaluations,  and  differen- 
tiations on  expressions  represented  as  SLIP  lists. 3,4,5  It  allows  the 
performance  of  such  operations  as  multiplying  the  expression  (1  + x)  by 
the  expression  (1  - x)  to  obtain  the  expression  (1  - x2)  rather  than 
some  numeric  value.  SYMBOLANG  is  one  of  several  formula  manipulation5 
or  algebraic  manipulation  systems  in  existence. 


SYMBOLANG  has  been  used  extensively.  It  has  been  employed  in 
generating  coefficients  for  Taylor  series,  coefficients  for  the  solution 
of  partial  differential  equations  with  boundary  values,  expansion  of 
determinants,  and  differentiation  of  complicated  expressions.  In  addi- 
tion to  the  many  subroutines  provided  in  SYMBOLANG,  the  user  has  the 
ability  to  define  new  functions  at  will  and  incorporate  them  into  the 
system. 


SYMBOLANG  is  an  extension  of  SLIP,  which  SYMBOLANG  uses  as  a list 
processing  system.  The  virtue  of  using  a list  processing  system  for 
embedding  a symbol  manipulator  is  that  the  space  needed  to  treat 


^apidus,  A.,  and  Goldstein,  M.  Some  Experiments  in  Algebraic  Manipula- 
tion by  Computer.  Communications  of  the  Association  for  Computing 
Machinery,  1965,  _8  501-508. 


2Findler,  N.  V.,  Pfaitz,  J.  L. , 
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New  York,  Spartan  Books,  1972, 


and  Bernstein , H.  J,, 
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305-387. 


Four  High-Level 
SYMBOL ANG. 


3Weizenbaum,  J. , Symmetric  List  Processor.  Communications  of  the 
Association  for  Computing  Machinery,  1963,  6_,  524-544. 

4Hirsc’hberg,  M.  A.  SLIP  for  the  BRLESC  II  Computer.  BRL  R 1731,  July 
1974 


5Findler,  N.  Y. , Pfaitz,  J.  L.,  and  Bernstein,  H.  J. , Four  High- 
Level  Extensions  of  Fortran  IV:  SLIP,  AMPPL-II,  TREETRAN,  SYMBOLANG. 

New  York,  Spartan  Books,  1972,'  1-82. 

5Sammet,  J.  E. , Survey  of  Formula  Manipulation.  Communications  of  the 
Association  for  Computing  Machinery,  1966, _9,  555-569. 
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expressions  expands  and  shrinks  over  a large  range.  A list  processing 
language  allows  the  allocation  of  space  needed;  unneeded  space  can  be 
returned  to  the  system.  In  SYMBOLANG,  it  is  the  user's  responsibility 
to  return  unneeded  space  to  the  system. 

SYMBOLANG  consists  of  seventy  two  subroutines  and  is  operational 
on  the  BRLESC  II  computer.  The  SYMBOLANG-SLIP  system  is  stored  on  the 
disc  and  is  accessed  by  use  of  the  following  statement: 

4-  t r mnn  P\/r  itrr’  hit 

Lupori  lc  uisl,  oilainu,  All 

SLIP  is  stored  on  the  disc  and  is  accessed  by  use  of  the  following 
statement: 


* rnMDTrp  nmr  cttd  a t t 

VjUI  U 1.  UL,  | UlUVj  } OUl  f r\ULj 

SLIP  is  written  primarily  in  FORTRAN  and  allows  the  programmer  to 
create  lists,  insert  and  delete  items,  and  scan  items  on  a list.  In 
SYMBOLANG  a SLIP  list  may  be  thought  of  as  parenthesized  group  of 
quantities.  Since  a SLIP  list  may  be  placed  on  a SLIP  list,  one  must 
take  care  never  to  place  a list  on  itself. 

We  will  assume  the  reader  has  a knowledge  of  FORTRAN  and  a working 
knowledge  of  SLIP  in  the  sequel  of  the  text. 

Expressions  are  represented  as  SLIP  lists  in  the  following  manner 
(using  the  Baekus-Naur  form): 

00 

< expression  >-*-({<  term  > } q) 

oo 

< term  > ■*-  (quantity  { < factor  > < power  > } r>) 


< factor  >-<-{<  expression  > } 


00 

0 symbol 


< power  > •*-  quantity  | < expression  > , 

where  the  left-hand  side  of  each  line  represents  the  class  of  objects 
being  defined;  the  angles  < > enclose  class  names;  braces  { } indicate 
that  the  enclosed  objects  are  to  be  included  in  the  definition  for  each 
number  of  repetitions  indicated  by  the  numbers  on  the  right  brace;  | 
separates  alternative  definitions;  quantity  is  the  class  of  floating 
point  quantities  other  than  zero;  and  symbol  is  the  class  of  expression 
symbols. 

In  keeping  with  Bernstein's  notation,  double  brackets  will  be  used 
to  denote  real  world  expressions,  and  commas  will  be  used  between  list 
elements. 
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is  empty. 


The  expression  [ [x] ] is  a nonconstant  expression  consisting  of  one 
term  with  coefficient  1,  a single  factor  1HX,  and  the  power  1.  It  is 
represented  as  the  list  ((1.,  1HX,  1.)). 
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In  general  a Hollerith  constant  of  from  one  to  ten  nonblank  char- 
acters is  used  to  represent  each  symbol  in  an  expression.  All  function 
arguments  are  represented  as  lists.  An  expression  is  broken  up  into  a 
sum  of  terms,  where  each  term  is  broken  up  into  a product  of  factors 
raised  to  powers.  The  trivial  coefficient  and  power  1 axe  also  intro- 
duced wherever  necessary.  Parentheses  without  a function  name  have  been 
treated  as  the  appearance  of  a special  function  name  that  has  arbitrarily 
been  translated  into  the  expression  symbol  3H1.*. 

B.  Ordering  and  Simplification 

In  SYMBOLANG,  identical  terms  differing  only  in  their  coefficients 
are  combined,  and  a 0 coefficient  causes  a term  to  be  dropped  from  an 
expression.  Identical  factors  differing  only  in  their  powers  are 
combined,  and  a 0 power  causes  a factor  to  be  dropped.  If  the  constant 
expressions  [[0]]  and  [[1]]  are  raised  to  any  power,  they  are  unchanged. 


SYMBOLANG  imposes  an  ordering  on  all  expressions  and  their  com- 
ponents. All  operations  preserve  this  ordering,  and  some  impose  it. 

The  ordering  imposed  is  built  up  from  orderings  on  numbers  and  expressio 
symbols.  Numbers  are  ordered  arithmetically  and  expression  symbols  lexi 
cographical ly  using  the  following  ordered  alphabet: 


0123456789ABCDEFG 
H IJKLMNOPQRSTUVWX 
YZ  + -*  / ()$=’*  [ ] : + < > 

with  blank  coming  before  all  other  characters.  Thus,  3H1.*  is  less  than 
1HA,  which  is  less  than  2 HA A which  is  less  than  1H$, 


Factors  within  terms  and  terms  within  expressions  are  reordered  to 
achieve  expressions  of  the  lowest  possible  ordering.  Terms  and  factors 
are  combined  and  dropped  where  possible  as  previously  defined. 


a hand  Croat sd  expression  may  be  simplified  by  calling  subroutine 
SMPL  (one  of  the  72  SYMBOLANG  subroutines). 


SMPL  (LIST)  reorganizes  the  expression  LIST  into  simplified  form 
returning  LIST  as  its  value.  The  operation  is  timewise  expensive  and 
should  be  used  sparingly. 
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The  list  of  available  space  must  be  initialized.  For  this  we  can 

use: 


DIMENSION  WORK  (5000) 

/"i  att  T\1  T*r»P  /"  r\  »/■  r nnnN 

CALL,  UN  1 1 AO  (nUK\,  OUUUj 
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up  a pointer  and  100  public  lists  in 


COMMON  AVSL,  X(100) 

X(98),  X (99) , and  X(100)  are  used  extensively  in  SYMBOLANG  and 
should  not  be  altered  by  the  user. 
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XSQ , and  XTQ  so  the  user  should  not  define  routines  or  labeled  common 


blocks  with  these  characters. 


In  addition  to  the  SYMBOLANG  subroutines,  the  system  requires  many 
of  the  SLIP  subroutines.  For  completeness,  the  author  has  made  avail- 
able the  entire  SLIP  package  for  use  with  the  system. 

Following  is  a list  of  all  the  entry  points  in  the  SYMBOLANG- SLIP 
package , 


ADD 

ADN  LWR 

ADVDL 

ADVDR 

ADVLEL 

ADVLER 

ADVLL 

ADVLNL 

ADVLNR 

ADVLR 

ADVLWL 

ADVLWR 

ADVSEL 

ADVSER 

ADVSL 

ADVSNL 

ADVSNR 

A r\Y  rf  r» 

nuvoK 

ADVSWL 

a nwoum 
rtuvoim 

ALOAD 

ALOG 

ATAN 

ROT 

UW  A 

BREAK 

CONI 

COS 

CP 

CPYTRM 

DELETE 

DERROR 

DLSRDR 

DVSUM 

DVTRM 

EQUAL 

E VALUE 

EXP 

GETCOE 

HITENT 

ICHR 

ID 

INHALT 

INITAS 

INITRD 

INLIST 

INLSTL 

INLSTR 

INSBST 

INSUBT 

INTARG 

INTENT 

INTGER 

IRADLS 

IRALST 

T TV  A TV  TV  TV 
1 

ITSVAL 

KOKE 

LANORM 

LCNTR 

LDATVL 

LDIF 

LEX ICQ 

LIST 

LISTAV 

LISTMT 

LISTON 

LNKL 

LNKR 

LOC 

LOCARG 

LOCT 

LOFRDR 

LOOKUP 

LPNTR 

LPURGE 

LRDRCP 

LRDROV 

LSCMPR 

LSPNTR 

LSQADD 

LSQAPL 

LSQCDR 

LSQCMP 

LSQCNM 

LSQCPY 

LSQCXP 

LSQDEF 

LSQDES 

LSQDSF 

LSQERR 

r rnp  a tv 

t cr\p\i  r 
r 

T COfTU 
ljovyvj  i i»i 

I COT  HD 
buyi  L/I\ 

T C O T M T 

i_iovy  ini 

T qot  rn 

LSQMEX 

LSQMNl 

f ^0MN2 

LSQMN3 

LSQMN4 

LSQMN5 

LSQMN7 

LSQMNL 

LSQMTM 

LSQNAM 

LSQ0U1 

LSQ0U2 

LSQ0U3 

LSQ0U4 

LSQ0U5 

LSQ0U6 

LSQ0U7 

LSQ0U8 

LSQOUT 

LSQPNT 

LSQRAZ 

LSQSBS 

LSQTRC 

LSQTYP 

LSQUMN 

LSQUNl 

LSQUN2 

LSQUN3 

LSQUN4 

LSQUN5 

LSQUN6 

LSQUN7 

LSQVAL 

LSQWL 

icc  rnv 
Ljour l 

j CTcn r 

Lt«J  1 byii 

r cr\mv 

1 1T1I\1\ 

T CTDDfl 

aa  i i i\w 

r \r r D\n 

1-1  V 1_1  1\  V X 

LVLRVT 

MADATR 

MADLFT 

MADNBT 

MADNTP 

MADOV 

MADRGT 

MAKEDL 

MAN  1 

MAN2 

MAN3 

MAN4 

MAN5 

MAN6 

MAN7 

g 


MANY 

MONO 

MRICLSS 

MRKLST 

MTDLST 

MTLIST 

NAMEDL 

NAMTST 

NtivBui 

NEiVTOP 

M T!7  T»F1  I A T 

MA  A 'T’T  fT 

XT T IAr  T T 

XT  T TT  CTT 

XUIT  PTn 

INLi  il  VftL 

IN  un  1 V L, 

IN  ULiC  Li  Li 

IN  U LiO  1 Ll 

INULiO  1 I\ 

NUMPY 

NXTLFT 

NXTRGT 

PARMT 

PARMT2 

POPBOT 

POPMID 

POPTOP 

P OUTER 

PRESRV 

PRIPUT 

PRLSTS 

PUTLST 

RCELL 

RDLSTA 

REED 

REEL 

RESTOR 

RITEA 

RITEF 

RITEI 

RITEIT 

RITEO 

SAME 

SBST 

SCHATL 

SCHATR 

SEQLL 

SEQLR 

SEQRDR 

SEQSL 

SEQSR 

SfclUlK 

SETIND 

SETRAY' 

CIJTM 
Ol  LX1N 

C TNT 
O 11N 

CUT)  T 
orir  Li 

cm  wc 

J>JL  V L. 

CATM 

OL^XIN 

snni  it 

SREED 

SRTRM 

.^TRITT  R 

STRIND 

SUB 

SUBSBT 

SUBST 

SUBSTP 

SUBT 

SUMPY 

TAN 

TANH 

TERM 

TOP 

TRCAL 

ZERO 

TRUNC 

TSTCON 

VALARG 

VISIT 

The  pattern  of  code  in  SYMBOLANG  is  usually  to  create  operands,  per- 
form expression  operations,  and  to  erase  unneeded  operands. 
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routines,  and  using  subroutines  as  if  they  were  functions,  In  order  to 
call  a FORTRAN  function  on  the  BRLESC  II  computer,  one  must  provide  an 
extra  argument  for  the  value  of  the  function.  This  argument  does  not 
need  to  be  accessed  by  the  user,  just  provided  by  him.  For  instance,  the 
function  XFUN(A,  B)  may  be  called  as  follows: 


CALL  XFUN  (C,  D,  TEMP) 

TEMP  would  contain  the  value  of  the  function. 


In  using  a subroutine  as  a function,  the  subroutine  must  have  at 
least  one  argument.  The  value  of  the  subroutine  cannot,  however,  be 
assumed  to  be  available  for  use. 


Functions  snd  subroutines  nre  often  deeply  nested  in  SYMBOLANG, 


II.  SYMBOLANG  SUBROUTINES 


A.  Basic  Set 


is  arbitrary,  but  it 
or  complexity. 


QVMRm  iwn  wi  i i Ko 

4 riuv  uru  ivj  nx  i i i/v 

overcomes  artificial 
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PTouoines  according  to  function 


(1)  FUNCTION  ADD  (LIST1,  LIST2) 


ADD  is  one  of  the  older  SY'MBOLANG  functions  whose  use  should  be 
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(2)  SUBROUTINE  ALOAD  (ARRAY.  VAR1 . VAR2 , 
K = 1,55. 


. . , VARK) , where 


ALOAD  stores  the  machine  address  of  the  variables  VARi  into  the 
array  elements  ARRAY(I).  A call  to  ALOAD  should  be  made  to  store  ad- 
dresses of  arguments  of  any  routine  that  uses  INTARG,  LOCARG,  or  VALARG. 
ALOAD  has  a variable  length  calling  sequence,  ALOAD  returns  the  argument 
count  in  the  event  an  argument  is  a floating  point  999  (used  where 
multiple  entries  are  not  provided).  The  argument  count  is  common  block 
MNYARG. 


(3)  BREAK  (LIST1,  SYMBOL,  LIST2,  LIST3) 

BREAK  is  one  of  the  older  SYMBOLANG  functions  whose  use  should  be 
avoided.  It  has  a newer  counterpart  LSOTRC  which  should  be  used  instead. 
BREAK  adds  to  the  expression  LIST2  those  terms  of  the  expression  LIST1 
that  contain  the  expression  symbol  SYMBOL  and  replaces  LIST2  with  this 
result.  The  remaining  terms  of  LIST1  are  added  to  the  expression  LISTS 
(which  is  changed  to  this  new  value). 


(4)  CPYTRM  (LIST1 , LIST2) 


CPYTRM  is  one  of  the  older  SYMBOLANG  functions  whose  use  should  be 
avoided.  CPYTRM  removes  the  top  element  from  list  or  expression  LIST1 
and  places  it  on  the  bottom  of  list  or  expression  LIST2.  If  LIST1  is 
empty,  no  change  is  made  and  CPYTRM  returns  0;  otherwise,  CPYTRM  returns 
1.  Correct  ordering  is  not  necessarily  preserved,  and  changes  to  LIST1 
may  affect  expressions  of  which  it  is  a subexpression. 

(5)  DVSUM  (LIST!,  SYMBOL,  LIST2) 


DVSUM  is  one  of  the  older  SYMBOLANG  functions  whose  use  should  be 
avoided.  It  adds  to  the  expression  LIST2  the  result  of  dividing  ex- 
pression LISTI  by  expression  symbol  SYMBOL.  LIST2  is  replaced  with  this 
result  which  is  also  returned  as  the  value  of  the  function. 


(6)  DVTRM  (LISTI,  SYMBOL) 


DVTRM  divides  expression  LISTI  by  expression  symbol  SYMBOL. 

(7)  EVALUE  (LISTI,  SYMBOL  1,  QUANT  1,  . . .,  SYMBOL  6,  QUANT  6) 


avoided. 


one  of  the  older  SYMBOLANG  functions  whose  use  should  be 
A number  of  other  routines  are  available  for  evaluation. 


EVALUE  evaluates  the  expression  LISTI  with  the  expression  symbols  SYMBOL 
i set  to  the  floating  point  constants  QUANTi . The  evaluate  returned  as 
the  value  of  the  function  also  replaces  the  expression  LISTI.  EVALUE 


■Pi  vp/1 


'TT £ 1 -t 
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iiAwu  calling  58C|uGTiC6  on  BRLESC  II* 

value  of  999  should  be  used  for  SYMBOL  i to  signify  that  it  and  the 
arguments  which  follow  it  are  not  to  be  considered  in  the  evaluations. 
Arguments  following  the  999  flag  must  be  set  to  zero  to  insure  proper 
calling  sequence  termination. 


(8)  GETCOE  (SYMBOL,  QUANT,  LIST1,  LIST2) 

GETCOE  is  one  of  the  older  SYMBOLANG  functions  whose  use  should  be 

3 VH  i A T t flHHc  t-ho  c*  vnro  ec  i rvrj  T T <o  -P-Pi  /-  i on  t-  -P  tV\£ 
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symbol  SYMBOL  raised  to  the  floating  point  constant  power  QUANT  in  the 
expression  LIST1,  LIST2  is  replaced  with  this  expression  which  is  also 
returned  as  the  value  of  the  function. 

(9)  HITENT  (LIST1 , SYMBOL,  LIST2 , LIST3) 

IlITENT  also  containing  entry  INTENT  is  one  of  the  older  SYMBOLANG 

fJITlf't  1 nrVir\ca  neo  chnill  H K a oma  -I  At*  A T *•  O +■ « « •**>■*  e*  -HAcnlf  A A 4 *«.  /-r 
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a level  1 substitution  of  the  expression  LIST1  for  the  expression  symbol 
SYMBOL  in  the  expression  LIST2  to  the  expression  LIST3.  INTENT  acts 

4 X 

as  does  HITENT  except  that  the  substitution  is  performed  on  all  levels. 
LSQSBS  is  a replacement  for  HITENT  in  the  new  code. 


(iO)  ICHR  (I ALPHA) 


ICHR  returns  a number  for  a letter  of  the  alphabet  or  the  digits 
0-9.  One  is  returned  for  A,  2 for  B,  . . . , 26  for  Z,  27  for  0.  . . .. 
and  36  for  9. 


(11)  INLIST  (VAR,  QUANT 1 , QUANT2) 


INLIST  returns  a list  lepreseuting  the  expression  found  on  logical 
unit  or  list  QUANT1.  The  variable  or  array  element  VAR  will  also  contain 
the  newly  created  expression,  INLIST  is  a means  of  translating  FORTRAN - 
like  expressions  into  their  internal  representation  on  lists.  Input  is 
handled  line  by  line  with  columns  7 through  72  scanned  for  an  expression. 
If  a $ is  encountered,  or  the  word  NOMOR  is  encountered  in  columns  1-5, 
or  if  columns  1-5  are  nonblank,  the  scan  ends.  Expressions  are  formed 
from  symbols  and  numbers  combined  by  the  arithmetic  operators  +,  -,  *, 

/,  and  **,  balanced  parentheses,  and  commas. 


On  occasion,  it  is  desirable  to  have  the  capability  of  inputting 
an  equation.  In  this  instance,  if  QUANT2  is  set  to  3HVAL,  it  is  possible 
to  input  an  equation,  with  the  "="  being  added  to  the  operator  list. 

When  this  option  is  desired,  a double  $ is  needed  to  terminate  the  input 
expression.  In  all  cases  where  an  "=”  appears  and  this  option  is  not 
desired,  QUANT2  should  be  set  to  the  floating  value  of  999.  In  addition, 
the  expressions  before  the  last  equal  sign  (in  case  there  are  more  than  1) 
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are  effectively  comments.  For  instance,  the  input  line  of  the  form 
E = M*C**2  — 10$  is  translated  into  the  valid  representation  of  [[10jj. 
Even  though  these  expressions  before  the  last  equal  sign  are  effectively 
comments,  they  must  conform  to  the  syntax  for  input  expressions, 

(12)  INSBST  (LIST1,  SYMBOL,  LIST2 , LIST3) 

INSBST  is  one  of  the  older  SYMBOLANG  functions  whose  use  should  be 
avoided.  It  returns  the  result  of  adding  an  all-levels  substitution  of 
the  expression  LIST1  for  the  expression  symbol  SYMBOL  in  the  expression 
LIST2  to  the  expression  LIST3,  LSOSBS  is  a better  choice  for  this  type 
of  substitution. 


(13)  INSUBT  (LIST1 , SYMBOL,  LIST2 , LIST3,  SYMBOL 1 , QUANT1 , 
SYMBOLS,  QUANTS) 


INSUBT  is  a function  of  a fixed  number  of  arguments.  The  argument 


1 i cf 


; m i 1 A 


f 1 no t in  r 


j}  a iiuoLuig  point  value  of  999  for  the  first 
symbol  SYMBOL!  not  used,  followed  by  zeros  for  all  remaining  arguments. 
INSUBT  returns  the  result  of  adding  an  all  levels  substitution  of  the 
expression  LIST1  for  the  expression  symbol  SYMBOL  in  the  expression  LIST2 
to  the  expression  LIST3  and  then  truncating  on  the  expression  symbols 
SYMBOL!  to  the  powers  given  by  the  floating-point  quantities  QUANTi.  The 
expression  LIST3  is  replaced  by  the  result.  INSUBT  is  one  of  the  older 
SYMBOLANG  routines  whose  use  should  be  avoided.  LSQSBS  and  LSQTRC  should 
be  used  to  perform  the  operations  described. 


(14)  INTARG  (ARGS,  INT) 


INTARG  returns  the  value  of  the  INTth  argument  ARGS  of  the  routine 
that  called  INTARG.  A call  of  ALOAD  must  be  made  before  INTARG  may  be 
used. 


(15)  KOKE  (ALPHA) 

KOKE  stores  ALPHA  on  an  internal  list  and  returns  0 if  it  has  not 
encountered  the  last  item  to  be  stored  and  1 if  it  has.  KOKE  is  called 
by  PUTLST  an  older  SYMBOLANG  function  whose  use  should  be  avoided. 

(16)  LDIF  (LIST,  SYMBOL) 

LDIF  returns  the  derivative  of  expression  LIST  with  respect  to  the 
expression  symbol  SYMBOL.  Only  explicit  functional  dependence  is 
considered. 


(17)  LEXICO  (LIST) 

LEXICQ  is  a lexicographical  sort  of  the  expression  LIST. 
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LISTON  is  a two  argument  call  to  INLIST  for  input  on  the  standard 
input  unit  (at  the  BRL  unit  S). 

(19)  LOCARG  (ARGS,  INT) 


LOCARG  returns  the  machine  address  of  the  INTt’n  word  of  ARGS. 
call  of  ALOAD  must  be  made  before  LOCARG  may  be  used. 


A 


(20)  LOOKUP  (ITYP,  SYMBOL,  LIST,  INT) 

LOOKUP  defines  the  expression  LIST  to  be  the  partial  derivative  of 
the  function  represented  by  the  expression  symbol  SYMBOL,  with  respect 
to  the  argument  INT,  an  integer  quantity.  ITYP  is  used  when  the  function 

If  ITYP  is  not  zero,  a user-defined 


is  determined  by  user  definitions. 

funrti 


(21)  LSCMPR  (LA,  LB) 


LSCMPR  compares  the  terms  LA  and  LB  with  their  coefficients  removed. 
Zero  is  returned  if  LA  and  LB  are  equal,  -1  if  LA  is  less  than  LB,  and 
+1  if  LA  is  greater  than  LB. 

(22)  LSQADD  (LIST1,  LIST2) 

LSQADD  returns  the  sum  of  expressions  LIST1  and  LIST2. 

(2  3)  LSQAPL  (NAM,  LIST,  QUANT) 


LSQArL  applies  the  function  NAM  to  expression  LIST  and  all  sub- 
expressions of  LIST  if  QUANT  is  any  value  except  floating  point  999. 


(24)  LSQCDR  (LIST,  SEQ) 


LSQCDR  substitutes  a user-defined  argument  in  expression  LIST  using 
sequence  reader  SEQ. 


ras'i  r<;nrMP  ri  a tr  ttydpi 
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LSQCMP  is  a lexicographical  comparison  of  LA  and  LB.  If  ITYPE  is 
3HEXP  then  LA  and  LB  are  expressions,  if  ITYPE  is  3HTER,  they  are  terms, 
and  if  ITYPE  is  3HFAC,  they  are  factors.  If  ITYPE  is  anything  else, 

LA  and  LB  are  terms  to  be  compared  as  if  their  coefficients  were  1.0. 
Note  that  the  sequence  readers  are  advanced  if  ITYPE  is  3HFAC.  Zero  is 
returned  if  LA  and  LB  are  identical,  -1  if  LA  is  lexicographically  less 
than  LB,  and  +1  if  LA  is  lexicographically  greater  than  LB. 


(26)  LSQCNM  (SYMB0L1 , SYMB0L2) 
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identical.  -1  if  SYMB0L1  is  lexicographically  less  than  SYMBOL  2 and 
+1  if  SYMBOL  1 is  lexicographically  greater  than  SYMB0L2. 


(27)  LSQCPY  (SEQ) 


LSQCPY 

r>  i i I'TUrtrt  4- 
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returns  a list  which  is  a copy  of  the  items  to  the  right  of 
setting  of  sequence  reader  SEQ. 


(28)  LSOCXP  (LIST) 


LSQCXP  returns  the  list  or  expression  LIST  or  a full  copy  of  LIST 
via  LSSCPY.  The  choice  is  determined  by  the  setting  of  the  variable 
LSQCSX  in  common  block  LSQCSX.  If  the  variable  is  0,  the  list  itself  is 
returned;  otherwise,  LSSCPY  (LIST)  is  returned. 

(29)  LSQDEF  (SYMBOL,  LIST1,  LIST2 , NOCUR) 

LSQDEF  defines  the  symbol  SYMBOL  to  be  the  expression  LIST1  in  which 
the  expression  symbols  on  the  list  LIST2  hold  the  places  into  which  any 
arguments  applied  to  SYMBOL  are  to  be  substituted.  Both  LIST1  and  LIST2 
are  erased  by  the  call.  The  user  really  uses  LSQDEF  to  define  his 
functions  which  may  be  evaluated  later  in  his  program. 


When  NOCUR  is  5HNOCUR,  the  expression  symbol  SYMBOL  is  made  to  appear 
undefined.  This  permits  SYMBOL  to  appear  within  LIST1  without  causing 
infinite  loops. 

LSQDEF  must  be  called  with  four  arguments;  however,  the  first 
argument  which  is  an  integer  999  effectively  terminates  the  calling 
sequence.  The  result  of  the  calling  sequence  in  this  instance  should  be 
padded  out  with  zeros. 


(30)  LSQDES  (LIST) 

LSQDES  erases  the  list  or  expression  LIST.  When  a list  is  no  longer 
needed,  it  should  be  destroyed  to  make  cells  available  for  other  lists. 

(31)  I .SODS  F (FUNCT.  SYMROLl 

V. J •% v- 7 J 


LSQDSF  aids  in  defining  and  evaluating  functions.  The  FORTRAN 
function  FUNCT  (defined  in  an  external  statement)  is  applied  to  the  next 
expression,  if  any,  reached  by  a right  advance  of  LSQDMF(13).  If  this 
expression  is  constant,  a constant  is  returned;  otherwise,  the 
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found  by  the  right  advance  of  LSQDMF(13)  as  its  argument. 


(32)  LSQERR  (ARG) 
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LSQERR  writes  an  error  message  and  terminates  the  SYMBOLANG  nan. 

ARG  is  a Hollerith  word  giving  the  name  of  the  routine  in  which  the  error 
occurred. 


f n>  *r  -w 
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LO^UAK  (CJXMBOL) 


LSQGAR  returns  the  next  expression  to  be  reached  by  a right  advance 

an  error  is  given  with  expression 


of  LS0DMF(13)  if  there  is  one.  If  not 
symbol  SYMBOL  as  the  error  message.  If  SYMBOL 
if  there  is  no  next  expression  in  the  advance. 


0,  a zero  is  returned 


(34)  LSQGNF  (SEQ,  NAME,  NOARGS) 

SEQ  is  a sequence  reader  assumed  pointing  to  the  first  item  of  a 

term.  SEQ  is  advanced  until  it  points  to  a non  sublist 
and  the  last  value  of  the  advancement  flag  is  returned  as  the  value  of 
LSQGNF.  NAME  and  NOARGS  contain  the  name  of  the  function  in  the  factor 
and  its  number  of  arguments. 

(35)  LSQGTM  (LA,  VAR,  POW1 , P0W2,  LVAL) 

See  LSQTRC  with  a calling  sequence  of  LA,  VAR,  POWl,  P0W2 , LVAL  +1, 

1. 

(36)  LSQIDR 

LSQIDR  is  a subroutine  called  by  LDIF  (on  its  first  use)  to  define 
the  derivatives  of  the  expression  symbols:  3H1.*,  3HSIN,  3HC0S , 3HTAN , 

3HC0T,  3HSEC , 3HCSC,  6HARCSIN,  6 HARCCSC,  4HSINH,  4HC0SH,  4HTANH,  4HCOTH, 
4HSECH,  4HCSCH,  7HARCSINH,  7HARCCOSH,  7HARCTANH,  7HARCC0TH , 7HARCSECH, 
7HARCCSCH,  3HEXP,  and  3HL0G.  The  defined  derivatives  are  respectively: 

1.,  COS(X).  -SIN(X),  SEC (X) **2 , -CSC(X)**2,  TAN (X) *SEC(X) , -COT  (X) *CSC (X) 
1/(1  - x**2)  **.  5 - 1/  (1-X**2)  **.  5 , 1/(1  + X**2)  , -1/(1  + X**2)  , 
X**-l*(x**2-l)**-.5,  -X**-l* (X**  2-1)  **-. 5 , .COSH(X) , SINH(X) , SECH(X)**2, 
-CSCH (X) **2 , -SECH(X) *TANH(X) , -CSCH(X) *C0TH(X) , 1/(X**2  + 1)**.5, 

1/ (X**2-l) **. 5,  1/ (1-X**2) , - 1/ (X**2-l) , -X**-l*(l-X**2) **-. 5, 
-X**-i*(X**2+I) **-.5,  EXP (X) , 1/X. 


l-5'  j 
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LSQINI  is  a subroutine  called  by  LSOVAL  (on  its  first  use)  to  define 
the  expression  symbols:  3HSIN , 3HCOS , 3HEXP,  3HL0G,  3HTAN , 6HARCTAN , 

4HTANH,  3H1.*,  2HV. , 2HQ.  and  3HIF.  The  definitions  represent 
respectively:  SIN(X),  COS(X),  EXP(X),  LOG(X),  TAN (X) , ARCTAN (X)  , TANH(X) , 

the  evaluates  of  the  arguments  enclosed  in  parentheses  (3H1.*),  the  argu- 
ment evaluated  twice  and  returned  as  its  value  (2HV. ),  the  unevaluated 
argument  (2HQ. ),  and  a conditional  expression  such  that  scanning  from 
left  to  right,  it  returns  the  first  even  numbered  argument  following  an 
odd  numbered  argument  that  evaluates  to  0 (3HIF.).  The  value  returned 
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is  an  unevaluated  argument.  If  there  is  an  odd  number  of  arguments, 
the  last  argument  is  returned  unevaluated.  In  all  other  cases,  0 is 
returned. 

Note  LSQINI  does  not  define  all  of  the  functions  whose  derivatives 
are  defined  in  LSQIDR. 

(38)  LSQLCO  (LIST) 

LSQLCO  returns  a lexicographically  ordered  copy  of  expression  LIST 

(39)  LSQMEX  (LIST1,  LIST2) 

LSQMEX  returns  the  product  of  expressions  LIST1  and  LIST2. 

(40)  LSQMNL  (OUANT 1 , . . , . OUANT25) 

LSQMNL  returns  a list  on  which  its  arguments  QUANTi,  i = 1,  25, 

appear.  LSQMNL  has  entries  L5QMN1,  LSQWJ2,  LSQMN3,  LSQMN4 , LSQNN5, 

and  LSQMN7  to  account  for  calls  with  1,  2,  3,  4,  5,  and  7 arguments. 

(41)  L-SQMTM  (LIST1,  LIST2) 

LSQMTM  returns  the  product  of  terms  LIST1  and  LIST2, 

(42)  LSQNAM  (LNAM) 

LSQNAM  returns  a unique  name  for  each  time  it  is  called  in  its 
argument  and  for  its  value. 

(43)  LSQOUT  (QUANTI,  ....  0UANT25) 

LSQOUT  outputs  items  for  function  LSQPNT.  The  argument  list  is 
indexed  through  and  obeyed  sequentially.  LSQOUT  has  QUANTi,  i = 1,  25 
arguments.  Entries  LSQ0U1,  LSQ0U2,  LSQ0U3,  LSQ0U4,  LSQOU5,  LSQOU6, 

LSQ0U7,  and  LSQ0U8,  are  provided  to  account  for  calls  with  1,  2,  3,  4, 

5,  6,  7,  and  8 arguments  respectively.  The  following  argument  blocks 
are  defined: 

INITIAL  ARGUMENT  NUMBER  OF  ARGUMENTS  IN  BLOCK 

5HFLUSH  1 

Forces  current  line  out. 

4HEDGE  2 

Sets  right  margin 

6HMARGIN  2 

Sets  left  margin  to  integer  represented  by  second  argument  of 
the  block. 


4 HUN IT 


2 

Switches  output  to  a new  logical  unit. 

5HCFLAG  2 or  3 

If  second  argument  of  the  block  is  zero,  no  special  continuation 
character  is  to  be  used  (i.e.,  in  initial  state);  otherwise,  the 
second  argument  specifies  a column  to  continue  a continuation 
flag  (usually  6)  and  the  left  most  non-blank  character  of  the 
third  argument  is  used  in  that  column  for  continuation  images. 

0 2 

Outputs  characters  to  a block  or  the  end  of  a word. 

Positive  integer  2 

Outputs  the  specified  number  of  characters  from  the  second 
argument  which  may  spread  over  several  machine  words. 

(44)  LSQPNT  (LIST,  LTITLE,  VALU) 

LSQPNT  is  used  to  print  expressions  LIST  and  titles  for  expressions 
LTITLE.  All  other  uses  of  LSQPNT  are  for  internal  use  of  SYMBOLANG. 

The  third  argument  should  be  set  to  floating  point  999. 

(45)  LSQRAZ  (LIST1,  LIST2) 

LSQRAZ  returns  the  result  of  raising  the  expression  LIST1  to  the 
expression  LIST2  power. 

(46)  LSQSBS  (LIST1 , SYMBOL,  LIST2 , INT) 

LSQSBS  returns  the  result  of  substituting  expression  LIST1  for  the 
expression  symbol  SYMBOL  in  the  expression  LIST2  through  the  level 
specified  by  the  integer  quantity  INT.  If  INT  is  -1,  a substitution  is 
made  on  all  levels  (an  infinite  loop  may  occur  if  LIST1  contains  SYMBOL). 

(47)  LSQTRC  (LIST,  SYMBOL,  QUANT1,  QUANT2 , INTI,  INT2) 

LSQTRC  returns  an  expression  derived  from  expression  LIST1  by  re- 
taining certain  terms  of  the  expression  through  level  INTI,  an  integer 
quantity,  and  all  terms  on  deeper  levels.  A term  is  retained  if  it 
contains  the  expression  symbol  SYMBOL  to  a nonconstant  power  or  to  a 
power  between  the  floating  point  quantities  QUANT1  and  QUANT2.  INT2 
helps  to  determine  how  the  truncation  occurs.  If  INT2  is  1 or  2,  one 
gets  SYMBOL**EXPRESSION  TRUNCATED.  If  INT2  is  2 or  4 one  gets  SYMBOL**P 
truncated,  P between  QUANT1  and  QUANT2,  and  for  INT  of  1 or  3,  one  gets 
SYMBOL**P  truncated  with  P not  between  QUANT1  and  QUANT2. 

(48)  LSQTYP  (QUANT,  VAR) 
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LSQTYP  returns  -1  if  the  quantity  QUANT  is  not  a list  and  0 or 
greater  if  it  is  a list  or  expression.  The  variable  VAR  is  set  to  the 
quantity  QUANT  if  QUANT  is  not  a list.  If  QUANT  is  to  be  considered  a 
list  but  not  necessarily  an  expression,  then  VAR  will  not  contain  useful 
information. 
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LSQUM1  removes  k (k  = 1,  25)  items  from  the  bottom  of  list  LIST, 
the  bottom  item  going  to  QUANTi,  then  next  to  QUANTi- 1,  etc.  The 
respective  entries  LSQUN1,'  LSQUN2,  LSQUN3,  LSQUN4,  LSQUN5,  LSQUN6,  and 
LSQUN7  are  provided  to  handle  calls  with  1,  2,  3,  4,  5,  6,  or  7 arguments. 
It  undoes  the  effect  of  MANY. 
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LSQVAL  returns  the  result  of  evaluating  expression  LIST  through 
level  INT,  an  integer  quantity.  If  INT  is  negative,  no  limit  is  in 
effect.  If  INT  is  0,  a copy  of  the  expression  LIST  is  returned.  The 
one  argument  call  (i.e.,  a call  with  INT  = 999)  is  equivalent  to  a call 
with  negative  INT. 


(51)  LSQWL  (LIST) 

LSQVVL  returns  a location  to  VISIT  in  order  to  evaluate  the  express- 
ion LIST.  The  scanning  of  arguments  is  made  possible  using  LSQDMF (13) 
as  a communication  cell. 

(52)  MANY  (LIST,  QUANTI,  ....  QUANT25) 

MANY  places  QUANTi,  i ~ 1,  25,  in  turn  to  the  bottom  of  LIST  and 
returns  LIST  as  its  value.  The  entries  MAN1,  MAN2,  MAN3,  MAN4,  MAN5, 
MAN6,  and  MAN7  are  provided  for  call  with  1,  2,  3,  4,  5,  6,  and  7 
arguments. 


(53)  MONO  (QUANT) 

MONO  provides  "the  quantity  QUANT  with  an  integer  name s 


(54)  NUMPY  (LIST,  QUANT) 


NUMPY  is  one  of  the  older  SYMBOLANG  routines  whose  use  should  be 
avoided.  It  replaces  the  expression  LIST  with  the  result  of  multiply- 
ing this  expression  with  the  floating  point  constant  QUANT  and  returns 
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(55)  POPMID  (VAR) 

POPMID  returns  the  list  element  on  top  of  which  the  variable  VAR 
is  sitting,  removes  that  element  from  the  list,  and  advances  VAR  to 
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the  left , VAR  is  treated  as  a sequence  reader  but 
it  was  newly  created. 


(56)  POWER  (LIST,  SYMBOL) 


POWER  returns  the  leftmost  power  to  which  the  expression  symbol 
SYMBOL  is  raised  in  the  expression  LIST,  considering  only  the  first 
level.  Zero  is  returned  if  SYMBOL  is  not  found. 


(57)  PRIPUT  (HOL,  LIST) 


PRIPITT  is  one  of  the  older  SYMBOLANG  routines  whose  use  should  be 
avoided.  PRIPUT  is  equivalent  to  LSQPNT  and  is  used  for  outputting. 


CCQl 


U 1 LjO  1 


( f T CT  AMAMTM 

i I , 


At*  ill  rr*A  n \ 

QUAiN  IZ  / ) 


PUTLST  is  one  of  the  older  SYMBOLANG  routines  whose  use  should 
be  avoided.  PUTLST  returns  the  expression  LIST  onto  which  the  terms 
QUANTi  have  been  added.  If  PUTLST  is  called  as  above  but  with  a leading 
0 in  the  calling  sequence,  PUTLST  will  not  reorder  factors. 

(59)  SAME  (QUANT) 


SAME  provides  the  quantity  QUANT  with  a floating-point  name. 

(60)  SBST  (LIST1 , SYMBOL,  LIST2 , LIST3) 


SBST  acts  as  does  INSBST  except  that  substitution  is  limited  to 
one  level. 


(61)  SETRAY  (A,  N,  QUANT) 

SETRAY  stores  the  quantity  QUANT  into  the  array  A whose  dimension 

is  N. 


(62)  SMrL  (LIST) 


SMPL  reorganizes  the  expression  LIST  into  simplified  form,  returning 
LIST  as  its  value.  The  reorganization  is  an  actual  change  and  can 
properly  alter  an  existing  expression.  The  call  is  a costly  one  whose 
use  should  be  made  sparingly. 


(63)  SOLVE  (LIST1 , SYMBOL,  LIST2) 


the  older  SYMBOLANG  routines  whose  use  should  be 
avoided.  It  assures  that  the  expression  LIST1  is  linear  in  the 
expression  symbol  SYMBOL  and  solves  the  equation  (LIST1  = 9)  for  SYMBOL. 
The  solution  is  added  to  expression  LIST  2 which  is  replaced  by  the  sum. 
The  sum  is  returned  as  the  value  of  the  function. 
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(64]  SRTRM(L] 

SRTRM  replaces  term  L with  term  L sorted  lexicographically. 


QI1R  ( T TQT1  T T CTO  ^ 

i.  y L1U14  ) 


SUB  is  one  of  the  older  SYMBOLANG  routines  whose  use  should  be 
avoided.  It  returns  the  result  of  subtracting  expression  LIST1  from 
expression  LIST2.  LIST1  is  destroyed  and  LIST2  is  replaced  with  the 
difference. 


(66)  SUBT(LIST,  SYMBOL,  LIST2 , LIST3,  SYMB0L1 , QUANT 1 , . . ., 
<;YMRnr ; 0UANT31 

SUBT  acts  as  does  INSUBT  except  that  the  substitution  is  limited  to 
one  level.  If  SYMBOLi  is  set  to  the  floating  point  value  of  999,  it 
ends  the  effective  calling  sequence  which  should,  however,  be  padded 
out  with  zeros. 


r a ~r\ 
La ' y 


T TCT1 

Ijl  U 1 ^ J 


QUANT5) 


SYMBOLI,  QUANT 1, 


cvunni  c 


o i 1‘iuulj 


SUMPY  is  one  of  the  older  SYMBOLANG  routines  whose  use  should  be 
avoided.  It  returns  the  result  of  adding  the  product  of  expressions 
L1ST1  and  LIST2  to  expression  LIST3  and  truncating  on  the  expression 
symbols  SYMBOL  i to  the  powers  QUANTi.  The  first  SYMBOLi  set  to 
floating  point  999  terminates  the  effective  calling  sequence;  however, 
the  calling  sequence  should  be  padded  out  with  zeros. 

(68)  TRCAL  (LIST,  SYMBOLI,  SYMB0L2) 

TRCAL  is  one  of  the  older  SYMBOLANG  routines  whose  use  should  be 
avoided.  It  finds  the  lowest  possible  power  to  which  the  symbol  SYMB0L2 
is  raised  in  those  terms  of  the  expression  LIST  that  contain  the 
expression  symbol  SYMBOLI.  Then,  the  expression  LIST  is  replaced  by  a 
trmicated  copy  in  which  those  terms  that  contain  the  expression  symbol 
SYMBOLI  to  a power  greater  than  that  minimal  power  of  SYMBOL  2 are 
discarded. 


(69)  TRUNC  (LIST,  SYMBOL,  QUANT) 


i cvunni  a Kin  i n aq  uU  acp  iic 

^ A 1 I'JUVJjTUlU  iUutlUtl)  nilUJV 


avoided.  It  replaces  the  expression  LIST  with  the  result  of  discarding 
all  first  level  terms  that  contain  the  expression  symbol  SYMBOL  to  either 
a power  greater  than  the  floating  point  quantity  QUANT  or  a non-constant 
power. 


TSTCON  returns  0 if  the  expression  LIST  is  not  a constant  expression. 
If  it  is  a constant,  it  returns  i and  the  variable  VAR  is  set  to  the 
equivalent  FORTRAN  constant  value. 

(71)  VALARG  (ARGS,  INT) 

VALARG  is  the  same  as  INTARG  except  that  it  has  a floating  point 

name. 


(72)  ZERO  (1,  J) 


a i . x ay 


B.  Expanded  Capability  Set 

+ ■ ■ - - 

In  addition  to  the  basic  SYMBOLANG  subroutines,  several  other 
routines  are  discussed  by  Dr.  Bernstein.  These  routines  are  presented 
below. 


ni  r cnmi  caiiamti 

LSOCON  returns  an  expression  equivalent  to  the  constant  QUANT. 

(2)  LSQGCO  (SYMBOL,  QUANT,  LIST) 

LSQGCO  returns  the  coefficient  of  the  expression  symbol  SYMBOL 
raised  to  the  power  QUANT  in  the  polynominai  LIST. 

(3)  SUBROUTINE  LSQIEQ 

LSQIEQ  causes  6HEQUAL  to  be  defined  as  a means  of  creating  defini- 
tions. 6HEQUAL  arises  with  the  three  argument  call  of  INLIST,  where 
one  wants  to  input  an  equation.  LSQIEQ  limits  6HEQUAL  to  two  useful 
arguments . 


f a ^ cuonniiTTMc  T QOTDM 

J JUDI\UUliiHC  Lov^iriN 

LSOIPN  defines  the  expression  symbol  5HPRINT  to  Drint  the  value  of 

■n.  i.  / x 

its  first  argument  with  its  second  argument  as  a name. 

(5)  LSQNDR  (LIST,  SYMBOL,  N) 

LSQNDR  returns  the  Nth  derivative  of  expression  LIST  with  respect 

+ rt  n i An  r*  < rmU  a 1 C Vll P O T 

L.U  CAplCOOXUU  Z>ytUU\J±  O I i'iUU  l_i  a 

(6)  LSONEB  (LIST,  N) 

LSQNEB  returns  the  Nth  element  from  the  bottom  of  list  LIST. 
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(7)  LSONET  ("LIST,  Ml 

LSQNET  returns  the  Nth  element  from  the  top  of  list  LIST. 

(8)  LSQXPN  (LIST) 


LSQXPN  returns  a version  of  the  expression  LIST  in  which  parentheses 
have  been  largely  removed.  Parentheses  remain  only  for  expressions  to 
positive  or  negative  fractional  powers  and  to  the  power  -1  unless  a non- 
constant power  is  encountered. 


III.  PROGRAM  EXAMPLES 


A.  Input  and  Output 


The  following  code  may  be  used  to  input  and  output  expressions: 

PROGRAM  INPUT 
DIMENSION  SP  (500) 

CALL  INITAS  (SP.  500) 

1 CONTINUE 

CALL  INLIST  (LA,  5HINPUT,  999,  TEMP) 

CALL  LSQPNT  (LA,  2HLA,  999. , TEMP) 

GO  TO  1 
END 

The  following  input  cards  run  with  the  above  program  produce  the 
output  shown  below. 

Input 

A1 1$ 

An  j.  a i i 

ruL  » rv  a J. 


NOMOR 


Output 

Aii$ 

T A — ATT 
Li /A  - All 

$ 

$ END  OF  EXPRESSION 

A12  + A1 1$ 

LA  = All  + A12 

$ 

/h  rim  Ar  r1  vnnrCr  Tf\M 

$ cunu  ur  nArKhboiuiN 
NOMOR 
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T A - A 
ijr\  -*  u 


$ END  OF  EXPRESSION 


An  equation  may  be  processed  by  changing  the  call  to  INLIST  above 
as  follows: 


CALL  INLIST 

The  following 
below: 


(LA,  5HINPUT,  3HVAL,  TEMP) 

input  combined  with  this  chsnge  produces  the  output 


$ 

fh 


E = M*C**2$$ 

LA  = EQUAL.  (E,  C**2*M) 


rvrn  nr 


cinu  ur 


rvnnrppT  AX  i 


EArKnOOlWIN 


Note  that  the  $■$  is  needed  to  properly  terminate  a legitimate  3 
argument  call  of  INLIST.  Note  also  the  use  of  the  cell  TEMP  with  a 
call  of  a FORTRAN  function. 


B.  Arithmetic  Calculation 


The  following  program  may  be  used  to  add  and  multiply  expressions: 
PROGRAM  TEST 


DIMENSION  SP 

(500) 

CALL 

INITAS 

(SP, 

500) 

CALL 

IN  LI ST 

(LA, 

5HINPUT,  999,  TEMP) 

CALL 

INLIST 

(LB, 

5HINPUT,  999,  TEMP) 

CALL 

T QODMT 

1 1 A 

rr  a 
v.  J-,/'  y 

2HLA, 

999.  , 

TEMP) 

CALL 

LSOPNT 

N 

(LB. 

2HLB . 

999.  , 

TEMP) 

LC  = 

LSQADD 

(LA, 

LB) 

LD  = 

LSQMEX 

(LA, 

LB) 

CALL 

LSQPNT 

(LC, 

2HLC, 

999.  , 

TEMP) 

CALL 

LSQPNT 

(LD, 

2HLD, 

999.  , 

TEMP) 

A A T T 

LALL 

rvTnn 

EA1  i 

FMn 

Lii  1 IV 


If  combined  with  the  input  below,  the  ensuing  output  results: 

INPUT 

All$ 

A12$ 

tf  CMn  TMDHT 
q>  j_.i  i lj  ui 

OUTPUT 
A1 1$ 

A12$ 
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LA  = All 


$ END  OF  EXPRESSION 
LB  = A12 

$ 

$ END  OF  EXPRESSION 
LC  = All  + A12 

$ 

ckth  cvnncccTAM 

1-.1XL/  Ur  l_i  AT  1YL.  X V_/l  i 

LD  = All  * A 12 

$ 

$ END  OF  EXPRESSION 
C.  Differentiation 


The  following  piOgram  allows  one  to  differentiate  expressions  wit 
respect  to  the  expression  symbol  X. 


PROGRAM  TEST 
DIMENSION  SP  (500) 

CALL  INITAS  (SP , 500) 

1 CONTINUE 

CALL  INLIST  (LA,  5HINPUT,  999,  TEMP) 

LB  = LDIF(LA,  1HX) 

CALL  LSQDES (LA,  TEMP) 

CALL  LSQPNT(LB,  5HDERIV,  999,,  TEMP) 

CALL  LSQDES (LB,  TEMP) 

GO  TO  1 
END 

Note  the  use  of  LSQDES  to  destroy  unneeded  expressions  and  thus 


The  input  below  produces  the  output  which  follows. 


INPUT 

1$ 

X$ 

yv 

X**N$ 

1/X**N$ 

SIN (A*X) /COS (A*X) $ 
F(X,X,X)$ 

L0G(N*X)  - N*L0G (X) $ 


OUTPUT 
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1$ 

DERIV=0 

$ 

$ END  OF  EXPRESSION 
X$ 

DERIV= 1 

$ 

$ END  OF  EXPRESSION 
X**2$ 

DERIV=2*X 

$ 

$ END  OF  EXPRESSION 
X**N$ 

DERIV=N*X**(-I+N) 

<t 

-f 

$ END  OF  EXPRESSION 
1/X**N$ 

DERIV=-N*X** (- 1-N) 

$ 

$ END  OF  EXPRESSION 
S IN (A*X) / COS (A*X) $ 

DERIV=A  + A*COS (A*X)**(-2)*SIN (A*X) **2 

$ 

$ END  OF  EXPRESSION 
F(X,X,X)$ 

DE RI V= PARTI AL(N. 0,1)  + PARTIAL(N. 0,2)  + PARTIAL(N. 0,3) 

$ 

N.O  = F (X, X, X) 

$ 

$ END  OF  EXPRESSION 
LOG (N*X)  - N LOG(X) $ 

DERIV=-N*X**f-l)  + X**(-l) 

$ 

$ END  OF  EXPRESSION 
NOMOR 

nnnTu-n 

UL  I\JL  V ” U 

S 

$ END  OF  EXPRESSION 
D.  Truncation 

The  following  program  may  be  used  to  truncate  a programmed  expres- 


PROGRAM  TEST 
DIMENSION  SP (500) 

CALL  INITAS (SP,  500) 

LA  = LSQMN3(LSQMN1(3. ) , LSQM'JSfS.,  1HX,  11.), 
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1 


The  program  above  produces  the  following  output: 

LA  = 3 + 5*X**11  - 4*X**22*Y 

$ 

$ END  OF  EXPRESSION 
LAI  = 3 

$ 

$ END  OF  EXPRESSION 
LA2  = 3 + 5*X**11 

$ 

$ END  OF  EXPRESSION 

LA3  = 3 + 5*X**11  - 4*X**22*Y 

$ 

$ END  OF  EXPRESSION 
LA4  = 3 + 5*X**11 

$ 

$ END  OF  EXPRESSION 
E.  Evaluation 


The  following  program  redefines  a definition  and  simplifies  an 
expression : 


PROGRAM  TEST 
DIMENSION  SP(10Q0) 

CALL  INITAS  (SP,  1000) 

LA  = LIST  (9) 

CALL  LSQDES  (LSQVAL(LA,  999),  TEMP) 
CALL  LSQDES  (LA,  999) 

CALL  LSQDEF  (3H1.*,  LSQml  (LSQMJ3 


( i i rvi i & i £ £ 

t A . , J.  . 


1 ^ x 

a.  ) ) , 


LSQNMI  (10H$.l.  *$.$.$), 


CALL  LSQDEF  ( 3HC0S , LSQMN 1 ( LSQNf! 4(1.,  LSQMN 2 ( LSQMN  1(1.), 
1 LSQMN3 (- 1 . , 3HSIN,  2.)),  3H.1*,.5)),  0,  999  TEMP) 

10  CALL  INLIST  (LA,  5 H INPUT,  999,  TEMP) 

IF  (LI STMT (LA)  . EQ. 0)  CALL  EXIT 
CALL  LSQPNT  (LA,  5HINPUT,  999,  TEMP) 

LB  = LSQVAL  (LA,  999) 


CALL  LSQDES 
CALL  LSQPNT 
LC  = LSQVAL 
CALL  LSQDES 
CALL  LSQPNT 
CALL  LSQDES 
GO  TO  10 
END 


(LA,  TEMP) 

(LB,  5HVALUE  1,  999. 
(LB,  999) 

(LB,  TEMP) 

(LC,  6HVALUE2 , 999., 
(LC,  TEMP) 


, TEMP) 


TEMP) 


The  following  inputs  for  the  above  program  produce  the  output  shown 
below. 


INPUT 

COS**2  + SIN**2$ 

COS (X) **2  + SIN (X) **2$ 

COS (X) *SIN (X) $ 

OUTPUT 

COS**2  + SIN**2$ 

INPUT  = COS* *2  + SIN**2 

$ END  OF  EXPRESSION 

VALUE  1 = SIN**2  + (1  - SIN**2) 

$ 

$ END  OF  EXPRESSION 
VALUE 2 = 1 

$ 

$ END  OF  EXPRESSION 

COS (X) **2  + SIN (X) **2$ 

INPUT  = COS (X) **2  + SIN (X) **2 

$ 

$ END  OF  EXPRESSION 

VALUE  1 = SIN(X) **2  + (1  - SIN**2,  X) 

$ 

$ END  OF  EXPRESSION 
VALUE 2 = 1 

$ 

$ END  OF  EXPRESSION 
COS (X) *SIN (X) $ 

INPUT  = COS (X) *SIN (X) 

$ 

$ END  OF  EXPRESSION 

VALUE 1 = SIN(X) * ( 1 - SIN**2 , X)**5E-1 

$ 

$ END  OF  EXPRESSION 

VALUE2  = (1  - SIN(X) **2) **5E- 1*SIN (X) 

$ 

$ END  OF  EXPRESSION 
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IV.  CONCLUDING  REMARKS 


The  full  power  of  SYMBOLANG  becomes  apparent  after  one  begins  to 
use  it.  The  subprograms  of  SYMBOLANG  combined  with  those  of  SLIP  give 
the  user  a wide  variety  of  options  for  solving  a large  class  of  problems. 
In  addition,  one  has  at  his  disposal  all  of  the  FORTRAN  language. 

SYMBOLANG  may  be  expanded  by  the  user  to  include  functions  and  sub- 
routines of  a general  or  specific  nature.  One  such  set  of  possible  ex- 
pansions is  an  integration  package  (see  e.g.  Slagle,  J.  R. , A 
Heuristic  Program  that  Solves  Symbolic  Integration  Problems  in  Freshman 
Calculus,  Journal  of  the  Association  for  Computing  Machinery,  1963,  10, 
507-520).  The  user  also  has  at  his  disposal  through  SYMBOLANG  the 
capability  of  defining  his  own  functions  and  incorporating  them  into  the 
system. 
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